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The gN final state is isospin-selective and thus provides access to the spectrum of excited nucleons 
without being affected by excited A states. To this end, the world database on eta photoproduction 
off the proton up to a center-of-mass energy of i5 ~ 2.3 GeV is analyzed, including data on differential 
cross sections, and single and double polarization observables. The resonance spectrum and its 
properties are determined in a combined analysis of eta and pion photoproduction off the proton 
together with the reactions -kN —> -kN, rjN, KA and KT,. For the analysis, the so-called Jiilich 
coupled-channel framework is used, incorporating unitarity, analyticity, and effective three-body 
channels. Parameters tied to photoproduction and hadronic interactions are varied simultaneously. 

The influence of recent MAMI T and F asymmetry data on the eta photoproduction amplitude is 
discussed in detail. 

PACS numbers: 11.80.Gw, 13.60.Le, 13.75.Gx. 


I. INTRODUCTION 

The determination of the spectrum of excited baryons 
from experimental data is necessary to understand Quan¬ 
tum Chromodynamics at low and medium energies. In 
this non-perturbative regime, quark models [1, 2] and 
lattice calculations [3-5] predict more exited states than 
found so far in partial-wave analyses of experimental 
data. This dilemma is known as “missing resonance prob¬ 
lem” [6]. In the past, the dominant source of informa¬ 
tion on resonance properties was provided by elastic ttN 
scattering [7-9]. The analysis of inelastic reactions is, 
however, essential [10, 11] when aiming at a reliable ex¬ 
traction of the entire spectrum and, consequently, at an 
identification of missing states that might couple predom¬ 
inantly to channels other than ttN. 

Among the inelastic channels accessible in ttN scatter¬ 
ing, the rjN channel plays a crucial role. It couples ex¬ 
clusively to states with isospin 1 = 1/2 allowing for the 
extraction of N* states unaffected by contributions from 
A* states. Moreover, the rjN channel opens at relatively 
low energies, in a region which is populated by numer¬ 
ous nucleon resonances. For example, some states like 
the four-star Af(1535)l/2“ resonance are known to have 
a large rjN branching ratio. Less well-established reso¬ 
nances like the Af(1710)l/2+, whose parameters are only 
weakly constrained from elastic ttN scattering [9, 11-13], 
may show a noticeable signal in their r]N decay [10, 11]. 
A narrow structure discovered in eta photoproduction on 
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the neutron [14-16] at around i? = 1.68 GeV could also 
appear in eta production on the proton target. 

In this respect it is unfortunate that the database for 
the reaction ttN —>■ rjN is problematic over the whole en¬ 
ergy range. Its coverage in scattering angles and energies 
is too limited to perform a well-founded resonance anal¬ 
ysis. Furthermore, for energies around 100 MeV from 
the threshold and beyond, only some experiments were 
performed and many of those are known to suffer from 
systematic uncertainties [9, 10, 17]. 

An alternative experimental window has opened in re¬ 
cent years with high-precision measurements of cross sec¬ 
tions and polarization observables in eta photoproduc¬ 
tion at photon-beam facilities like ELSA, JLab or MAMI, 
see Refs. [18-20] for reviews. For example, the first mea¬ 
surements from the JLab FROST target (polarization 
E for 7 p —>■ TT+n) have appeared only recently [21], and 
many more polarization data are expected. The database 
for eta photoproduction is not yet as large as the one for 
pion photoproduction, but it is rapidly growing. In addi¬ 
tion, the data already available are of much higher qual¬ 
ity than those for the pion-induced reaction ttV — > rjN. 
Recently, the first data for the beam-target asymmetry 
F in 7 p —>■ Tjp was presented by the A2 collaboration at 
MAMI, together with a measurement of the target asym¬ 
metry T [22]. 

A key question is whether eta photoproduction data 
can be used to access the poorly known N* branching ra¬ 
tios to the TjN channel. The combined analysis of elastic 
ttN scattering and pion photoproduction determines, at 
least in principle, the helicity couplings and ttN branch¬ 
ing ratios. With known helicity couplings, then indeed 
the data on eta photoproduction allow one to pin down 
the resonance rjN branching ratios, without having to re- 
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sort to the problematic ttN —>■ r]N data [9, 10, 17]. This 
argumentation is, however, limited by the fact that even 
for pion photoproduction the database is not yet fully 
complete. The extension of the analysis to eta photo¬ 
production will then improve the knowledge of the rjN 
branching ratios, but, due to incomplete databases, will 
also lead to changes both in the helicity couplings and 
branching ratios, as will be seen. 

Given the discussed mismatch in data quality, it is 
not recommendable to only fit model parameters tied 
to photoproduction and leave the hadronic interaction 
unchanged: The hadronic amplitude, poorly fixed from 
the ttN —>■ r]N data, appears as a sub-process in photo¬ 
production. Thus, to avoid any bias, an unconstrained 
fit is needed. In this way, the higher statistical weight 
of the eta photoproduction data even provides a better 
constraint on the hadronic ttN —>■ rjN amplitude. 

The photoproduction of 77 mesons is also a prime can¬ 
didate for a “complete experiment”. From a mathemati¬ 
cal point of view, a complete experiment [23] consists of a 
set of eight carefully chosen observables, which resolve all 
discrete ambiguities up to an overall phase [24, 25]. For 
example, a complete set including F and T is given by 
{a, S, T, P, E, F, Cx, Ox} [24]. For experiments with real¬ 
istic uncertainties, however, eight observables are not suf¬ 
ficient [26-28]. Less than eight observables are required 
in a truncated partial-wave analysis [29, 30]. 

By contrast, for a complete experiment on the reac¬ 
tion Tr~p rjn only four observables are needed [31, 
32]. Given the data situation for that reaction, a re¬ 
measurement would greatly advance our understanding 
of the r]N final state. Physics opportunities with a pion 
beam are discussed in Ref. [33]. In any case, the current 
database does not contain a complete set of observables, 
neither for photon- nor pion-induced eta production. One 
then has to resort to other approaches that often combine 
data from different initial and final states. 

Over the years, a variety of theoretical approaches has 
been applied to analyze the pion- and photon-induced 
production of Jj mesons. For example, photoproduction 
of r] mesons in the resonance region not too far from 
threshold was studied in the framework of unitarized 
chiral perturbation theory in Refs. [34-38]. Consider¬ 
ing a broader energy range up to and beyond 2 GeV, 
iF-matrix [7, 39-44] and unitary isobar [45] models are 
practical tools to perform an analysis of large amounts 
of data. Sometimes, the real part of the self-energies 
is neglected and only on-shell intermediate states are 
maintained, which reduces the complexity of the calcu¬ 
lations. For the purpose of a combined analysis of dif¬ 
ferent reactions over a wide energy range, so-called dy¬ 
namical coupled-channel (DCC) models provide a partic¬ 
ularly suited framework. Theoretical constraints of the 
5'-matrix, like two- and three-body unitarity, analyticity, 
left-hand cuts and complex branch points, are manifestly 
implemented or at least approximated. This enables the 
reliable determination of the resonance spectrum in terms 
of pole positions, residues, and helicity couplings in the 


complex energy plane. The production of ry mesons in 
DGG approaches was studied, e.g., in Refs. [46, 47]. 

Here, we extend the Jiilich model, a DGG approach 
pursued over many years [10, 48-51] starting with 
Ref. [52] , to perform a simultaneous analysis of the pion- 
induced reactions ttN —>■ ttN, 'K~p —>■ rjn, 

, TT+p —?► and the photon-induced reactions 

7 P —>■ TT^p, TT+n, and pp. We allow the hadronic ampli¬ 
tudes themselves to vary, in addition to the parameters 
tied to photoproduction. As discussed, this is necessary 
because the quality of the data in •n~p —>■ pn is much in¬ 
ferior to the data in eta photoproduction. In a simultane¬ 
ous fit to all pion- and photon-induced data, we observe 
that, indeed, the eta photoproduction data have a strong 
influence. This influence reaches beyond the electromag¬ 
netic resonance properties and affects also resonance pole 
positions and hadronic branching ratios. 

The paper is organized as follows: In Sec. II, a short 
overview of the applied formalism is given. For a more 
detailed introduction of the semi-phenomenological ap¬ 
proach to meson photoproduction we refer the reader to 
Ref. [51]. In Ref. [10], an extensive description of the 
hadronic Jiilich DCC framework is provided. In Sec. Ill A 
we describe the data analysis and in Sec. Ill B the fit re¬ 
sults are shown. The extracted resonance parameters 
are presented and discussed in Sec. IV. Technical details 
about the renormalization of the nucleon mass are sum¬ 
marized in an appendix. 


II. FORMALISM 

In the approach (referred to as “Jiilich model”), the 
hadronic scattering potential is iterated in a Lippmann- 
Schwinger equation formulated in time-ordered pertur¬ 
bation theory (TOPT) and two-body unitarity is, thus, 
automatically fulfilled. The three-body tittN states are 
parameterized through the channels pN, aN and ttA. 
These effective three-body channels are included dynam¬ 
ically, i.e., the tttt and ttTV subsystems match the cor¬ 
responding phase shifts [49]. The analytic structure of 
the amplitude is given through real and complex branch 
points [ 12 ] and the real, dispersive contributions of the 
intermediate states. Moreover, t- and rt-channel ex¬ 
changes of known mesons and baryons constitute the 
non-resonant part of the amplitude and serve as “back¬ 
ground”. While the u-channel diagrams approximate the 
left-hand cuts, t-channel meson exchanges are essential to 
achieve three-body unitarity [53]. Note that the latter is, 
at the moment, only approximately satisfied in the Jiilich 
model. By means of this explicit treatment of the back¬ 
ground, strong correlations between the different partial 
waves and a non-trivial energy and angular dependence of 
the observables are generated. Although t- and w-channel 
processes are necessary for analytic structure and uni¬ 
tarity, they do not fully determine the amplitude. Bare 
resonance states are included as s-channel processes. In 
contrast to previous versions of the approach, here we 
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also allow for additional contact interactions. Such inter¬ 
actions do not spoil the analytic properties ensured by s-, 
t- and u-channel interactions. They absorb physics be¬ 
yond the explicit processes and, thus, increase the model- 
independence of the approach at the cost of a few more 
parameters. Practically, the changes in the amplitudes 
induced by the contact terms are comparatively small 
and the so-called background is still dominated by the t- 
and u-channel exchanges. Details are given in the follow¬ 
ing section. Note that contact terms are also included in 
Ref. [32]. 

In Ref. [50] the approach was extended to the 
strangeness sector incorporating the final state 

in the analysis. In the Julich2012 model of Ref. [10], 
the spectrum of nucleon and A resonances was extracted 
from a simultaneous analysis of the reactions ttN —>■ ttN, 
r]N, KA and KTi. The extension of the Jiilich ap¬ 
proach to pion photoproduction in a field-theoretical for¬ 
mulation, that respects the generalized off-shell Ward- 
Takahashi identity, was achieved in Ref. [54]. In Ref. [51], 
by contrast, the photon interaction is approximated in a 
phenomenological framework and the Julich2012 analysis 
serves as final-state interaction. The flexible formulation 
of Ref. [51], used to study the world data on pion photo¬ 
production on the proton, proved to be capable of analyz¬ 
ing large amounts of data while at the same time main¬ 
taining the analytic properties of the Jiilich approach. 
This framework, with the addition of contact terms, will 
be applied in the present study. 


A. Pion-induced reactions 

The pion-induced reactions are treated within the 
Jiilich dynamical coupled-channel formalism [10]. The T- 
matrix which describes the scattering process of a baryon 
and a meson can be formulated in the partial-wave basis 
and reads 


Tf,„iq,p',E) = Vf,^{q,p\E) 

OO 

+ X] / dpp^Vf,i^{q,p,E)G^{p,E)T^^{p,p',E) . (1) 


In Eq. (1) and in the following, E always means the 
scattering energy in the center-of-mass (c.m.) frame and 
q = ]g I {p' = \p'\) is the modulus of the outgoing (incom¬ 
ing) three-momentum that can be on- or off-shell. The 
channel indices v, p and k, represent the incoming, out¬ 
going and intermediate meson-baryon pairs, respectively. 

The propagators G„ for channels with stable particles, 
i.e. K = ttN, rjN, KA, or ATE, are given by 


Gk (p, E) e - Ea (p) - Eb (p) -I- ie 


( 2 ) 


with Ea = \J + p2 and Eb = being the on- 

mass-shell energies of the intermediate particles a and b 


in channel k with respective masses nia and mb- In case 
of the channels with unstable particles pN, aN and ttA 
that parameterize the ttttN channels, the propagators are 
more complex; for details see Ref. [49]. 

The scattering potentials V^i, can be decomposed into 
a pole and a non-pole part 


NP 


= k;. 


NP 


= V 

* 


E 

2=0 


. -A/ . 

Ifi',i ly,i 

E — rrfi 


( 3 ) 


The quantity denotes the sum of all t- and w-channel 
exchange diagrams, while comprises the s-channel 
resonance graphs. The functions correspond to 

the bare creation (annihilation) vertices of a resonance i 
with bare mass and are constructed from an effective 
Lagrangian which can be found in Table 8 of Ref. [50]. 
Explicit expressions of y)].j and y“.j are also given in 
Ref. [50], cf. also Appendix A of Ref. [10]. The ex¬ 
change potentials that constitute also derived from 
effective Lagrangians, are compiled in Appendix B of 
Ref. [10]. The decomposition of Eq. (3) is slightly modi¬ 
fied in the present approach. We implement here contact 
terms that do not introduce any singularities and that 
are used to absorb physics not explicitly contained in the 
parameterization through s-, t- and w-channel processes. 
The contact terms are introduced in a separable form and 
separately for every partial wave. 


^CT ^ CT;a CT;c 

uy I a ly 

^ rriN 


( 4 ) 


where the (y^"^’®) have the same functional form 

as the resonance vertices y)].j (y)).^) in Eq. (3). The asso¬ 
ciated couplings in the y^"*" are now new free parameters 
to be adapted in the fit to the data. 

Formally, the numerator structure of Eq. (4) is the 
same as the one of s-channel pole terms, to ensure the 
correct threshold behavior. Also, the contact terms carry 
channel indices. Formally, we can treat the contact terms 
as bare s-channel processes and absorb the contributions 
in the definition of V^: 


^y ' 


E 

2=0 


CT 


A/ • A^ - 

I Y 

E-my 


( 5 ) 


For a compact notation, we will no longer distinguish 
between bare resonance vertices y and structures of the 
contact term y®"'’ in the following. Then, the index i can 
refer to a bare s-channel resonance vertex y or one of the 
terms in the numerator of Eq. (4) . 

Similar to the potential the scattering matrix 
can also be written as the sum of a pole and a non-pole 
part, 

= + ( 6 ) 

with the unitary defined through 


^NP 
^ ^y 




nNP 


K, 


( 7 ) 
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Here and in the following, we do not display function 
arguments and integration symbols in favor of a more 
compact notation. Besides the resonance poles arising 
from s-channel diagrams in T^, a dynamical generation 
of poles in is also possible, as explained in Refs. [10, 
48]. 

Note that the pole part can be evaluated from 
the non-pole part . To this purpose one defines the 
dressed creation (annihilation) vertex r°.j (r“.j) via 


r° = 

Mi* 


■pa 


^ij 


0-° - 4- -7° ■ G 

y 

-7“ + yNP Q a 

!~ / ^ ly;i 5 

y 


( 8 ) 


where S is the self-energy. The indices i and j label the 
s-channel states or a contact diagram in a given partial 
wave. 

The pole part is then given by 

p _ pa n .. pc ('q'l 

^I 11 > ~ ^ fj,;i ^13 ^ v,J 

with the resonance propagator Dij. For example, if there 
are two s-channel resonances with masses mi and m 2 
(indices i, j G { 1 , 2 }) plus one contact term (indices i, j = 
3), the expression reads 


D-^ 





•pa 

/pa pa pa \ 

= r^;2 



Vr^;3/ 

'E-m\- 

Ell —S12 

-Sl 3 \ 

— E21 

E - m^- E22 
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— E31 

— E32 

miv — E33 / 


( 10 ) 


The decomposition into pole and non-pole part per¬ 
formed here has mostly technical reasons: The numer¬ 
ical evaluation of the non-pole part is much more time- 
consuming than the evaluation of the pole part. This 
leads to an effective, nested fitting workflow as discussed 
in detail in Refs. [10, 51]. Here, we have introduced the 
contact terms technically on the same footing as reso¬ 
nances, which allows to fit those terms computationally 
more effectively. 

In Sec. 2.2 of Ref. [10], the renormalization of the nu¬ 
cleon mass and coupling in the presence of two bare s- 
channel states in the Pn partial wave was derived. With 
the implementation of contact diagrams, this procedure 
has to be extended to two resonances and one contact 
diagram. This is addressed in Appendix A. 


B. Photon-induced reactions 

A field-theoretical description of the photoproduction 
amplitude within a gauge-invariant framework that re¬ 
spects the generalized off-shell Ward-Takahashi iden¬ 
tity [55-57] was successfully applied in the analysis of 


pion photoproduction in Ref. [54]. In this study an ear¬ 
lier version of the Jiilich model was utilized to provide 
the hadronic final-state interaction. This field-theoretical 
method allows one to gain insight into the microscopic 
reaction dynamics of the photo-interaction. By contrast, 
the semi-phenomenological approach to pseudoscalar me¬ 
son photoproduction, developed in Ref. [51] and used 
here, is more flexible and facilitates the analysis of a 
large amount of data. Here, the photo-interaction ker¬ 
nel is approximated by energy-dependent polynomials, 
while the hadronic final-state interaction is provided by 
the Jiilich DCC model described in Sec. HA. The for¬ 
malism is inspired by the GW-SAID CM12 parameter¬ 
ization [58] and will be applied in the present study. 
Nonetheless, we consider the present analysis as an inter¬ 
mediate step towards an expansion of the field-theoretical 
framework of Ref. [54]. A detailed introduction to the 
semi-phenomenological approach was given in Sec. 2.2 of 
Ref. [51]. In the following, we recapitulate the basic ele¬ 
ments. 

The multipole amplitude of the photoproduction pro¬ 
cess is given by 


+ J2Tf..{q,P,E)G,ip,E)V^^{p,E) , ( 11 ) 

where the index 7 denotes the initial jN state. is the 
hadronic half-off-shell T-matrix introduced in Sec. HA 
with the intermediate (final) meson-baryon channel k 
(/ i) and the corresponding off-shell momentum p (on- 
shell momentum q). Integration over the intermediate 
off-shell momentum p similar to Eq. (1) is suppressed 
here in the notation, following the convention of Eq. (7). 
In the present analysis the photon is allowed to couple 
to the intermediate channels k = ttN, iqN and ttA, while 
we have ttTV and pN as final states p. 

The photoproduction kernel Vfj,^ is written as 

V^^{p, E) = {p, E) + ^ -- . (12) 


Here, stands for the photon coupling to the non-pole 
part of the photoproduction kernel. The tree-level cou¬ 
pling of the jN channel to the nucleon and A resonances 
is represented by the vertex function 7 “.,, where i de¬ 
notes the resonance number in a given partial wave. The 
hadronic resonance annihilation vertex 7 “,^ is exactly the 
same as in Eq. (3), which results in the cancellation of the 
explicit singularity at E = m\. Both quantities, and 
7 l(.j, are approximated by energy-dependent polynomials 

P', 


AA^) 



(13) 


In Eq. (13), the vertex function 7 “ has the same form 
as 7 “ j in Eq. (12) but without any dependence on the 
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resonance number i. The polynomials P, for a given 
multipole, are parameterized as 


Fit A Fit B 


Pf{E) 

P^^{E) 


Y.9I 

^fj. 

JIC 

3=0 



'E-Eg 


V toat 

NP 

fP-Eg 

M.i 

V ITT-N 


^ g-Af 


(14) 


with g and A > 0 being free parameters that are fitted to 
the data. The upper limits of the summation indices j, 
ii and are chosen so as to permit a good data descrip¬ 
tion, but are restricted to be less than 4. In order to ful¬ 
fill the decoupling theorem, which states that resonance 
contributions are parametrically suppressed at threshold, 
the summation for P^ starts with j = 1. The expan¬ 
sion point Eg is chosen to be close to the nN threshold, 
Eg = 1077 MeV. In this way, the factor ab¬ 

sorbs the potentially strong energy dependence at the yfV 
threshold, which is not too far from the ttN threshold. 
Moreover, this factor guarantees a well-behaved multi¬ 
pole amplitude in the high-energy limit, although a more 
quantitative matching to Regge amplitudes remains to 
be done [59]. 

In order to achieve a good description of the high- 
precision data for pion photoproduction close to thresh¬ 
old, we take into account some isospin breaking effects, 
i.e. we apply different threshold energies for the tt'^p and 
the TT+n channels, as explained in Sec. 2.3 in Ref. [51]. 
Similarly, we take the physical threshold of the gp final 
state when calculating observables. Note that, in general, 
isospin-averaged masses are used in the Jiilich model. 

A multipole decomposition of the photoproduction am¬ 
plitude of pseudoscalar mesons can be found in Ap¬ 
pendix A of Ref. [51]. 


III. RESULTS 

A. Database and free parameters 

The database of the present study comprises the 
hadronic data used in Ref. [10]. This represents the world 
database on the reactions ttN gN, KA, and KY, up 
to E ^ 2.3 GeV, plus the nN —>• nN WI08 energy- 
dependent solution of the GWU/INS SAID group [60]. 
In addition, we include almost all published data on pion 
photoproduction off the proton up to £1 ^ 2.3 GeV, with 
some forward regions at high energy excluded [51]. Third, 
in this study we add the world database for the reaction 
7 P —>■ gp, again up to if ~ 2.3 GeV. These data were 
taken from the GW-SAID database [60]. 

In a first fit, called fit A in the following, we include 
all data except the recent eta photoproduction measure¬ 
ment of the transverse target asymmetry T and the 
beam-target asymmetry E by the A2 collaboration at 


ttA ^ ttA 
7r“p —^ gn 
n~p —>■ K^A 
n-p^ K°T,° 
n~p —> 

TT + p^ K+T,+ 


PWA GW-SAID WI08 [60] 
da/dfl, P 
da/dtl, P, P 
daldtd, P 
da j dPl 

dajdtd, P , P 
~ 6000 data points 


gp —^ tA'p 

da/dtl, E, P, 

T, Anai, G, H 

gp —>■ 7r’*’n 

dajdO,, E, P, 

T, Angi, G, H 

7P Tfp 

da/dtd, P, E 

dajdtl, P, E, T, F 


29,392 data points 

29,680 data points 


TABLE I: Data included in fits A and B. In contrast to fit 
A, fit B contains the recent MAMI measurements of T and 
F [22]. 


MAMI [22]. These data are added in a second fit, fit 
B. Performing these two fits allows us to estimate the 
influence of the new polarization data on the extracted 
resonance spectrum. An overview of the data included 
in the fits can be found in Table I. 

Gompared to the high-precision data nowadays avail¬ 
able in case of pseudoscalar meson photoproduction, the 
data situation for the pion-induced reactions is difficult in 
large parts, the lack of polarization measurements being 
one of the major issues. In Sec. 3 of Ref. [10] the situ¬ 
ation for the individual hadronic channels was discussed 
in detail. In the present study, we adopt the system¬ 
atic errors and mainly also the special weights applied to 
certain data sets. 

Moreover, we continue along the lines of Ref. [51] and 
apply an additional systematic error of 5 % to all photon- 
induced data in order to account for discrepancies in the 
data. The amount of 5% is an estimate at this point. 
More advanced techniques have been applied, e.g., by 
the GWU/SAID group allowing for normalization cor¬ 
rections [9, 58]. We plan to improve our analysis along 
these lines in the future. To compensate for the smaller 
number of data points for T and E in eta photoproduc¬ 
tion, those data are weighted in fit B with generic factors 
around 10 as found necessary for obtaining satisfactory 
fit results. To achieve a good description of the data 
at higher energies, additional weights have to be applied. 
The situation is similar in case of the beam asymmetry S 
in gp gp. Gompared to the number of data points for 
the differential cross section, 5680, only a few are avail¬ 
able for S, namely 189. Furthermore, we did not attempt 
to achieve a good description of the recoil polarization P 
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for 7 p —> r]p as only seven data points are available and 
their influence on the fit is very limited. The data are, 
however, included in the fit but no special weights were 
applied. 

In the Jiilich approach the free parameters tied to 
the hadronic interaction are the bare coupling constants 
and masses of the s-channel resonances, the strengths of 
the contact terms, and the cut-off parameters in the t- 
and u-channel diagrams constituting . In addition, 
there are certain couplings in some of the latter diagrams 
that cannot be connected to other coupling constants via 
SU(3) flavor symmetry and have to be fitted to data as 
well. In the present study we do not alter the pa¬ 
rameters but employ the values found in the Julich2012 
analysis of pion-induced reactions [ 10 ]. 

We do, however, vary the parameters tied to hadronic 
resonances and contact terms. This is a necessity in the 
current situation, given that the data of pion-induced 
eta production is of less quality than the eta photopro¬ 
duction data. The higher-quality photoproduction data 
can, thus, help to constrain the hadronic amplitude. In 
both fits A and B the number of s-channel resonance pa¬ 
rameters amounts to 128. For each of the 11 genuine 
7=1/2 and 10 genuine 7 = 3/2 resonances (i.e., those 
resonances included in form of a pole diagram in the po¬ 
tential [10, 50]) the parameters are given by one bare 
mass and the couplings to the channels ttN, pN, rjN, 
ttA, KA and TfS as allowed by isospin. Contact terms 
introduced in Eq. (4) are switched on in both fits in the 
S'!! and Pii partial waves for the ttN and pN channel, 
and in the P 13 partial wave for the ttN, pN, ttA, KA 
and KT, channel, giving rise to 9 additional fit param¬ 
eters. The free parameters that are used to tune the 
interaction of the photon with hadrons are the resonance 
parameters gfj and Af and the non-pole couplings 
and AJJ'J with p = ttA, ttA or TyfV, cf. Eq. (14). For¬ 
mally, all fit parameters up to j = 3 are implemented 
in the computer code. However, the actual number of fit 
parameters is chosen as required by data and the remain¬ 
ing parameters are set to zero. In order to achieve a good 
description of T and P in fit B, additional fit parameters 
that were set to zero in fit A had to be released. Thus, we 
have 443 parameters tied to the photon interaction in fit 
A and 456 in fit B. In total, the number of fit parameters 
adds up to 580 in fit A and 593 in fit B. 

The free parameters are adjusted to the data in simul¬ 
taneous fits of all pion- and photon-induced reactions 
using MINUIT on the JUROPA supercomputer at the 
Forschungszentrum Jiilich. 


B. Fit results 

In the following, only selected fit results for the reac¬ 
tions 7 p —>■ gp and 7 r“p gn are shown. Data sets 
with energies that differ by less than 5.5 MeV are some¬ 
times displayed in the same graph. The full fit results for 
all pion- and photon-induced reactions included in this 


analysis can be found online [61]. 

The definition of the various photoproduction observ¬ 
ables in terms of the multipole amplitudes is given 
in Appendix B of Ref. [51]. The convention agrees with 
the one of the SAID group [62, 63]. 


1. gp ^ gp 

In Figs. 1 and 2 we show selected results for the dif¬ 
ferential cross section. In the threshold region only very 
small differences can be observed between fits A and B. 
Starting at P ^ 1600 MeV, the two fits show some dis¬ 
crepancies at extreme angles, fit A providing a slightly 
better description of the data at medium energies, cf. 
E = 1629 MeV in Fig. 1. At higher energies, the differ¬ 
ences are most apparent at very forward angles. However, 
the data situation does not allow for an assessment re¬ 
garding which of the fits is best. This can be seen, e.g., at 
E = 2006 MeV in Fig. 2 where both fits seem to describe 
the data equally well. 

The situation is similar in case of the beam asymmetry 
S in Fig. 3. Disagreements of the fits A and B are visible 
predominantly at higher energies and at forward angles 
where there are no data. 

Only seven data points are available for the recoil po¬ 
larization P, cf. Fig. 4. Although the results of our two 
fits are very different, especially at higher energies, both 
describe the data more or less well but with some defi¬ 
ciencies. A larger database of this observable may help 
constrain the partial wave content. 

The fit results for the transverse target asymmetry T 
and the beam-target asymmetry P can be found in Fig. 5. 
These data were only included in fit B, meaning that the 
blue dashed line in Fig. 5 (fit A) represents a prediction 
for these observables. Note that older data for T from 
Ref. [77], that are partially in conflict with the MAMI 
data, were not fitted. The prediction of fit A at lower en¬ 
ergies E < 1600 MeV is acceptable for T and good for P, 
while the shortcomings of the prediction are considerable 
at higher energies. However, once the data are included 
(fit B) a good description over the whole energy range is 
achieved. 


2. n p ^ gn 

In Fig. 6 we show selected fit results for the differential 
cross section of the reaction tt~p — >■ gn. In addition to 
fits A and B, we display results from the old fit A of the 
Julich2012 analysis [10], in which only pion-induced re¬ 
actions were analyzed. We call the latter fit Ahad in the 
following, where the subscript serves as a reminder that 
the fit in Ref. [10] included only hadronic data. Since 
the resonance vertex function 7 “^^ appears in the con¬ 
struction of the hadronic amplitude in Eq. (3) and also 
in the photoproduction kernel in Eq. (12), the bare pa¬ 
rameters of 7 “jy have influence on the pion- as well as on 
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FIG. 1: Differential cross section of the reaction 7 p —>■ r)p. Dashed (blue) line: fit A; solid (red) line: fit B; data: MClO [64], 
BA07 [65], KR95 [66], NA06 ]67], DY95 [68], CR05 [69], DU02 [70], CR09 ]71], WI09 [72]. In this and all following figures, the 
numbers in the plots specify the pertinent center-of-mass energy i5[MeV]. 


the photon-induced production of the r]N final state. As 
much less data are available for tt~p —> rjn compared to 
7P —>■ r/p, the photon-induced eta production does impose 
constraints on the fit results of the hadronic r/N channel. 

As can be seen in Fig. 6, starting already at energies 
of if ~ 1.5 GeV, data sets from different experimental 
groups show conflicting behavior. Moreover, at energies 
E > 1800 MeV all data for the differential cross section of 
Tr~p — >■ rin stem from a measurement [78] deemed prob¬ 
lematic due to a miscalibration of the beam momentum, 
see Ref. [17] for details. Those data enter the fit with a 
much reduced weight and all three fits yield different re¬ 
sults. This is most apparent at higher energies, where fits 
A and B are dominated by photon-induced data. Here, 
fits A and B mostly differ at extreme forward angles while 
for other angles they coincide reasonably well. This dis¬ 
crepancy has been also noted for the differential cross 
section and beam asymmetry in fits A and B in eta pho¬ 


toproduction, mostly due to less precise data at forward 
angles. 

In Fig. 7 we present results of the recoil polarization 
in 7r“p —>■ r]n. In principle, the only published data for 
this observable [79] exhibit the same problems as the dif¬ 
ferential cross sections of Ref. [78] because the same ex¬ 
perimental set up was used. Those data are fitted with 
a very low weight, too. As in the case of the differential 
cross section, the differences in the three fits are obvious. 
However, in view of the discussed quality issues with the 
available data for tt~p —>■ rjn, we see no reason to enforce 
a better data description. 

The comparison made here demonstrates the need for 
a pion beam to re-measure the reaction 7r“p —> rjn. Note 
that only four observables are needed for a complete ex¬ 
periment, as discussed in the Introduction. 
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FIG. 2: Differential cross section of the reaction 7 p —>■ rjp. Dashed (blue) line: fit A; solid (red) line: fit B; data: cf. Fig. 1 and 
SU09 [73]. 
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FIG. 3: Beam asymmetry of the reaction 7 p —>■ rjp. Dashed (blue) line: fit A; solid (red) line: fit B; data: BA07 [65], ER07 [74]. 
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FIG. 4: Recoil polarization of the reaction 7 p — rjp. Dashed 
(blue) line: fit A; solid (red) line: fit B; data: HE70 [75], 
H071 [76]. 


C. Multipoles 

In Fig. 8 the electric and magnetic multipoles from fits 
A and B for the reaction 'yp —>■ rjp are shown. Addi¬ 
tionally, we display the multipoles of the Bonn-Gatchina 
BG2014-02 solution [85]. Note that those amplitudes 
were not included in the fit. 

Generally speaking, the multipoles extracted in our 
analysis and the ones of the Bonn-Gatchina group exhibit 
large differences. An exception is the i?o+ and, to a cer¬ 
tain degree, also the Mi_ multipole. While for dif¬ 
ferences between our fits A and B are hardly noticeable at 
all, the Mi_ multipole features clearly visible deviations 
in the two fits at energies higher than E ~I650 MeV. As 
will be discussed in Sec. IV A, this is the energy regime 
where the pole of the V(1710)l/2+ resonance is located. 
Among the lower multipoles, the Mi+ shows the most 
striking discrepancies between fits A and B, see also the 
discussion on the influence of variations in the P 13 par¬ 
tial wave on the description of T and E in Sec. IV A. 
For higher multipoles fit B sometimes shows a stronger 
energy dependence than fit A, cf. i? 4 _, M 4 - and E 4 +, 
M 4 +. In summary, the new MAMI data for T and F 
have a large impact on the multipole amplitudes. 

We observe that for lower partial waves, the eta pho¬ 
toproduction multipoles of the present study exhibit 
less agreement with the Bonn-Gatchina multipoles than 
our pion photoproduction multipoles with the Bonn- 
Gatchina BG2014-02 [85] or the GW-SAID CM12 solu¬ 
tion [58]. This suggests that the multipole content of the 
reaction yp —>■ pp is much less established than in the 
case of pion photoproduction where the various analyses 
agree better. Figures showing the pion photoproduction 
multipoles can be found online [61]. 


IV. RESONANCE SPECTRUM 

A resonance state is uniquely defined by its pole posi¬ 
tion in the complex energy plane, the residues associated 
with the channel transitions, and the Riemann sheet the 
pole is located on. With the exception of the physical 
sheet of the lowest lying channel, the poles can appear 
on various Riemann sheets, but not all of them are of 
physical interest. Usually, only the poles on the sheet 
which is closest to the physical axis are considered. We 
select this sheet by rotating the right-hand cuts of all 
channels in the direction of the negative imaginary E 
axis. In this way, we define the second sheet where all 
poles extracted in the present study lie. See Ref. [49] for 
a detailed discussion. 

In order to determine the pole positions, the scatter¬ 
ing amplitude has to be continued to the second Rie¬ 
mann sheet. For this purpose we apply the method of 
analytic continuation following Ref. [49] where the am¬ 
plitude on the second sheet is accessed via a contour de¬ 
formation of the momentum integration. The calculation 
of the residues proceeds via the formalism illustrated in 
Appendix G of Ref. [50]. Definitions of the normalized 
residue and the branching ratio into a specific channel 
are given in Sec. 4.1 of Ref. [10]. In case of the latter 
two quantities we use the same definitions as the Particle 
Data Group [ 86 ]. For a reliable extraction of the reso¬ 
nance parameters, the correct structure of branch points, 
including the complex branch points of the channels in¬ 
cluding unstable particles ttA, aN and pN, is crucial. 
In Ref. [12] it was shown that the absence of the latter 
might lead to false resonance signals. 

The definition of the photocouplings at the pole ^poig, 

= ^poiee*”' , (15) 

can be found in Appendix C of Ref. [51] and is identical 
to the definition given in Ref. [87]. The photocoupling 
at the pole characterizes the coupling of the pN channel 
to a resonance, independently of the final state in the 
reaction under consideration. Note that, in general, the 
complex Apgjg cannot be compared to the real-valued he- 
licity amplitudes , see Sec. D of Ref. [51] for further 
remarks. 

In Tables II to IV we list the pole positions, residues 
and the photocouplings at the pole of the present study. 
In addition to the values extracted in the current fits A 
and B we list the pole positions and residues found in fit 
A of the Julich2012 analysis [10], called fit Ahad in the 
present study, and the photocouplings of fit 2 from the 
Julich2013 analysis [51]. Note that in the latter study, 
the parameters of the hadronic T-matrix were not al¬ 
tered, i.e. the resonance pole positions and hadronic 
residues are the same as in fit Ahad of Ref. [10]. An 
overview of the pole positions of fit A and B is also given 
in Fig. 9. 

In Table III the ttA channel labeled ( 6 ) corresponds 
to the case where \ J — L\ = 1/2 and the one labeled (7) 
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FIG. 5: Transverse target asymmetry T and beam-target asymmetry F in the reaction 'yp pp. Dashed (blue) line: prediction 
of fit A; solid (red) line: fit B; data: Ref. [22]. 


to |J — L| = 3/2. Also the pN channel can couple to a 
resonance with a given in multiple ways, cf. Table 11 
in Ref. [10]. Here, we only quote normalized residues 
for ttA, since at energies well above the ttA threshold 
this channel can be regarded as being composed of the 
two stable particles tt and A. In general, the resonance 
coupling at the pole to a channel like ttA is a function 
of the center-of-mass momentum of the stable particle 
(that equals the summed momenta of the decay prod¬ 
ucts of the unstable particle). Here, we do not quote 
this function of gc.m. but choose gc.m. as on-shell three- 
momentum of a stable A of mass m = 1232 MeV and a 
pion. Obviously, this prescription does not lead to mean¬ 
ingful results for the very broad a in the aN channel, or 
the pN channel. In the latter channel, most resonances 
are not far above the threshold that is situated around 
E = (1.7 — i 0.075) GeV, and the p cannot be considered 
a stable particle. 


A. Discussion of specific resonances 

Compared to the earlier analysis of Ref. [10], the exten¬ 
sion of the model to eta photoproduction did not require 
the inclusion of additional bare s-channel states, and we 
find no new dynamically generated resonances either. In 
particular, there is no need to include a narrow state at 
around E k, 1.68 GeV. The narrow structure discovered 
in eta photoproduction on the neutron [14-16] is absent 
in the present analysis of eta photoproduction on the 
proton. 

In the following, when discussing selected resonance 
states, we always refer to the values quoted in Tables H 
to IV. 

Sii: While the real part of the pole position of the 
N(1535) l/2~ is very stable throughout all three fits 
A, B and Ahad, the width is by 30 MeV larger in the 
new fits that include eta photoproduction data. The new 
value is close to the one found in a recent analysis [88] of 
the GW-SAID WI08 solution [60], where elastic ttN and 
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FIG. 6: Differential cross section of the reaction n~p —^ rjn. 
Dashed (blue) line: fit A; solid (red) line: fit B; dash-dotted 
(green) line: fit Ahad of Ref. [10]. Data: filled circles from 
Ref. [80]; filled squares from Ref. [81]; empty triangles up from 
Ref. [82]; stars from Ref. [83]; empty squares from Ref. [84]; 
empty diamonds from Ref. [78]. Note that the data situation 
for this reaction is problematic, see text. 


ttN rjN data were fitted. Also the mass in the latter 
analysis is very similar to our fits. We obtain the same 
values for the normalized r]N residue in all three fits. 
While the elastic ttN residue is larger in the new fits, 
the coupling to this channel is still considerably smaller 
than the coupling to r]N. In both current fits A and B, 
the magnitude of the photocoupling is more than 

twice as large as in the previous fit 2 of Ref. [51] where 
only pion photoproduction data were considered. This 
change is related to the increase of the width for the 
A^(1535) 1/2+, because the resonance width and size of 
the photocoupling at the pole are strongly correlated. 
On the whole, comparing the resonance parameters of 
the iV(1535) 1/2“ of the older fits to the new ones, the 
inclusion of eta photoproduction seems to have noticeable 
impact for this resonance. The influence of the new T 
and F data [22], on the other hand, is rather limited 
as the parameters in fit A and B do not exhibit major 
differences. This observation is in agreement with the 
similarity of fit A and B for T and F at A ~ 1.5 GeV (cf. 
Fig. 5 ). 

By contrast, the pole of the second resonance in the 
Fii partial wave, N(1650)1/2- ■ , is located in an energy 
region where the deficiencies of the prediction of fit A for 
T and F become more apparent. Accordingly, slightly 
larger variations are found in the pole positions of fits 
A and B, and also compared to fit Ahad- Although the 
width is smaller in the new fits A and B, the resonance is 
still broader and has a higher mass than in Ref. [88]. In 
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FIG. 7: Polarization of the reaction 7r“p — >■ rjn. Dashed (blue) 
line: fit A; solid (red) line: fit B; dash-dotted (green) line: fit 
Ahad of Ref. [10]. Data: Ref. [79]. Note that the data situation 
for this reaction is problematic, see text. 


Ref. [89], however, the results from a Laurent-Pietarinen 
(L-l-P) expansion of the GW-SAID CM12 solution [58] 
for pion photoproduction give a pole position of Eq = 
1655(11) — i63.5(8.5) MeV, which is closer to our values. 
As in fit Ahad j the current fits reveal a strong coupling to 
the KY channels. The magnitude of the photocoupling 
to the Ai(1650) 1/2“ is more than twice as large in the 
fits including eta photoproduction compared to the older 
fit 2 of Ref. [51], where only pion photoproduction was 
considered. While fit 2 yielded a value smaller than the 
one of Ref. [89] , the photocoupling is now larger in fit A 
and B. 

In summary, the pole positions of the Af(1535) 1/2“ 
and the A^(1650) 1/2“ come closer to the values of the 
GW-SAID analysis. As resonance widths and sizes of 
helicity couplings are correlated, the values of the lat¬ 
ter also approach the values of that analysis. The eta 
photoproduction data have a strong influence on the Fn 
resonance properties. 

Pii: Besides the nucleon pole, we find two resonances 
in the Pn partial wave. One of them, the Roper reso¬ 
nance Af(1440) 1/2+ is dynamically generated from the 
interplay of the t- and rt-channel diagrams. As the fit pa¬ 
rameters corresponding to these diagrams are not 
altered in the present study the extracted resonance pa¬ 
rameters do not change much. The situation is differ¬ 
ent for the third state, the A^(1710) 1/2+. This explicit 
s-channel state was introduced in the Jiilich model in 
Ref. [10] mainly to improve the description of the pion- 
induced r]N and KK channels. Since it couples only 
weakly to the ttN channel its resonance parameters are 
poorly constrained from t:N elastic scattering. As can be 
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FIG. 8: Electric and magnetic multipoles for the reaction 7 p —>■ pp. (Black) dash-dotted line: BG 2014-02 solution [85]. Dashed 
(blue) line: fit A; solid (red) line: fit B. 


seen in Table II, the extension of the fit to new, inelastic 
reaction channels results in noticeable changes in the pole 
position. Moreover, also the inclusion of new observables 
for a specific reaction, here T and f in 7 p —>■ rjp in fit B, 
leads to significant variations not only in the pole posi¬ 
tion but also for the residues and photocouplings. The 
latter observation suggests that additional information 
from inelastic channels, e.g. in form of new polarization 
measurements, might help to fix the parameters of the 
1V(1710) 1/2+. 

In all our fits, the 7V(1710) 1/2+ has a lower mass and 
is narrower than in recent analyses by the ANL-Osaka 
{Eq = 1746 — il77 MeV) [46] and the Bonn-Gatchina 
groups [Eo = 1687 ± 17 - i(100 ± 12.5) MeV) [90]. In 
the L-l-P analysis of the GW-SAID CM12 solution in 
Ref. [89], a broad state with a higher mass associated 
with this resonance is found that can be alternatively ex¬ 
plained as the pN complex branch point. However, the 
authors state that additional information from other de¬ 
cay channels beside ttN is needed to distinguish between 
the two options. See also Ref. [ 88 ] by the same authors 


where the same conclusion was drawn. Note that in the 
present study the pN complex branch point is included 
explicitly. 

In addition to the V(1440) 1/2+ and the V(1710) 
1 / 2 + we find non-conclusive indications for another, very 
broad and dynamically generated pole at i? ~ 1.75 GeV. 

P 13 : We include one bare s-channel state in the P 13 
partial wave, the V(1720) 3/2+. Although we observed a 
noticeable sensitivity of the description of the rjN channel 
on variations in the Pl3 partial wave, the pole position 
of the V(1720) 3/2+ is very similar in fits A and B. The 
impact of the new T and F data from MAMI can be 
seen in the photocouplings of this state (fit B vs. fit A). 
This is reflected in the discrepancies observed in the Mi+ 
multipole in Fig. 8 . 

In different GW-SAID solutions [ 88 , 89] and in the 
Bonn-Gatchina analysis of Ref. [90] the iV(1720) 3/2+ 
has a pole position with a real part 20 to 80 MeV lower 
than in our fits and an imaginary part more than 50 MeV 
larger. By contrast, the ANL-Osaka group [46] finds val¬ 
ues closer to ours. 
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FIG. 9: Pole positions of the isospin /=l/2 (above) and 
/=3/2 (below) resonances extracted from fit A (filled squares) 
and from fit B (empty diamonds). For a better differentiation, 
the names and pole positions of some resonances are marked 
with different colors. The green crosses denote the branch 
points of the amplitude. Note that all cuts, starting at the 
branch points, are chosen in the negative Im E direction. 


We tested the influence of a second explicit resonance 
state in the P 13 partial wave but observed no significant 
improvement of the fit results. In this partial wave, the 
Bonn-Gatchina group finds strong evidence for a state 
named iV(1900) 3/2+ in the photoproduction of KK and 
KYi [90, 91]. It has also been confirmed in 7 p —>■ A'+A 
in an effective Lagrangian model [92] and in a covariant 
isobar-model single channel analysis [93]. It remains to 
be seen, whether this state is also needed in the Jiilich 
approach once the analysis is extended to kaon photo¬ 
production. Note that the A^(1900) 3/2+ is also in¬ 
cluded in the ANL-Osaka analysis [46] and in the Giefien 
model [94]. 

D 13 : While the real part of the pole positions of the 
A^(1520) 3/2“ is unchanged in fits A and B, the imagi¬ 
nary part is about 10 MeV smaller in fit B. In the previous 


fit Ahad the real and the imaginary part were similar to 
fit A, -2Im i?o=110 MeV. In the GWU-SAID solutions 
analyzed in Refs. [ 88 , 89] and the Bonn-Gatchina analy¬ 
sis [90], widths of about 110 MeV are also found. With 
-2Im Eo=78 MeV, a smaller width was extracted in the 
ANL-Osaka analysis [46]. 

Moderate changes in our three fits can also be observed 
for the values of the residues and photocouplings. Al¬ 
though the iV(1520) 3/2“ is well determined from elastic 
ttN scattering and no new information from this chan¬ 
nel was included in the new fits, certain changes in the 
resonance parameters are not surprising. Due to the well- 
known SD-wave interference in the pion-induced rjN pro¬ 
duction resulting in a rt-shape form of the differential 
cross section (cf. Fig. 6 ), the 7V(1520) 3/2“ shows some 
sensitivity to the parameterization of the rjN channel. As 
can be seen in Fig. 6 at i? = 1509 MeV and especially 
at E = 1576 MeV, the description of the data differs in 
all three fits. The energy bin at E = 1576 MeV is, on 
the other hand, a prime example for the systematic prob¬ 
lems in the data. Data at backward angles are available 
with small error bars, but not in agreement with other 
data spanning the entire angular region. Underestimated 
normalization problems can obviously change the angu¬ 
lar dependence significantly and have a large impact on 
the partial-wave content. Better data are called for. 

Di 5 , F 15 : Although the poles of the V(1675) 5/2“ 
and the V(1680) 5/2+ are located in an energy region 
where the prediction of fit A for T and E becomes worse 
(cf. Fig. 5), the pole positions and residues exhibit only 
minor differences in fit A and B. Still, as can be seen 
in Fig. 10, in the current fit B, the His and F 15 are 
important to achieve a good description of the new T 
data in eta photoproduction. Whereas for E a qualitative 
description of the data is feasible with the S'!!, P 13 , and 
ZI 13 partial waves alone, in case of T all S-, P-, D-, and 
F-waves are needed at medium and higher energies. 

In the current form of the approach, only one bare 
s-channel state is incorporated in the F 15 partial wave. 
For a discussion of a possible second explicit resonance 
we refer the reader to Sec. 4.3 of Ref. [10]. 

F 17 : The r]N normalized residue of the 7V(1990) 7/2+ 
is small. Still, the inclusion of the new polarization data 
for the reaction yp —>■ 7 ]p results in a pole position with 
the real part 53 MeV and the imaginary part almost 
60 MeV smaller in fit B than in fit A. The magnitudes of 
photocouplings, on the other hand, show much less vari¬ 
ations in the two new fits compared to the previous fit 2 . 
As remarked in Sec. 4.2 of Ref. [10], from an analysis of 
elastic ttN scattering not much evidence can be claimed 
for this resonance. However, in our current fit B the F 17 
partial wave seems to play a certain role in the param¬ 
eterization of T in 7 p —?> 77 P, cf. Fig. 10. At energies in 
the range of the pole position of the iV(1990) 7/2+, the 
Fi 5 alone plus the S-, P-, and D-waves does not yield 
a qualitative description of the data. However, evidence 
for this resonance from the current database is weak in 
general. 
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FIG. 10: Partial wave contribution to T (upper row) and F 
(lower row) in fit B. Solid (red) line: fit B; dashed (green) 
line: Sn, P 13 , 7 )i 3 wave only; dash-dotted (indigo) line: all 
S-, P-, D-, and F-waves. 


responding single-energy solutions do not provide enough 
information for correlated fits. Furthermore, also for 
eta photoproduction certain data sets were included with 
a higher weight. This renders an error estimation as per¬ 
formed in Ref. [51] impracticable for the present analysis. 

A comprehensive statistical error analysis is compli¬ 
cated by the large number of data points and free param¬ 
eters, typically inherent in the kind of analysis at hand. 
Such an analysis is, to our knowledge, not pursued in any 
of the current DCC approaches and we postpone a rigor¬ 
ous error analysis to future work. Without such an error 
analysis, the assessment of the significance of certain less 
well-determined states, like a potential A^(1750) 1/2+, is 
not possible. Concerning the significance of resonance 
signals, the systematic elimination of states as pursued 
in Ref. [9] is also a necessary task which we postpone to 
future work. 


V. SUMMARY 


Gi7, Gi 9, Hig; We include bare s-channel states 
identified with the A^(2190) 7/2“, Af(2250) 9/2“ and 
W(2220) 9/2+ resonance. The parameters of those broad 
states are less stable, as was already observed in the 
Julich2012 analysis [10]. 

Isospin 7 = 3/2 resonances: Since the rjN channel 
does not couple to resonances with I = 3/2, the pole 
positions and hadronic residues are very similar in fits 
A and B. The inclusion of eta photoproduction data can 
change the 7 = 3/2 resonances only indirectly through 
the mixed-isospin ttN channels in pion photoproduction, 
and the mixed-isospin channels in the reaction ttN —>■ 
TsTS. Nonetheless, the mass of the A(1930) 5/2“ is about 
40 MeV higher and the widths about 70 MeV larger in fit 
B. Smaller differences in the parameters from fits A and 
B can also be seen in the A(2200) 7/2“ and the A(2400) 
9/2“. The photocouplings at the pole are marginally less 
stable than the pole positions. 

Comparing fit A and fit Ahad of Ref. [10], the influence 
of pion photoproduction data included in the former fit 
but not in the latter is visible in the results for some 
states, as e.g. the A(1910) 1 / 2 +. Note that also the 
A(1232) 3/2+ changes its pole position slightly. 

In the analysis of pion photoproduction within the 
Jiilich framework [51] the uncertainties of the extracted 
photocouplings were estimated from re-fits based on dif¬ 
ferent re-weighted data sets. In Ref. [51], all data in¬ 
cluded in the fit entered with a universal weight of one. 
In the present study, however, the situation is different. 
As described in Sec. Ill A, the quality of the hadronic 
data requires a specific weighting of the various data sets. 
Moreover, in case of the elastic ttN channel we fit to the 
energy-dependent partial-wave amplitudes of the GW- 
SAID group for which no errors are provided. As a side 
remark, it should be noted that the error bars of the cor- 


Over the last years, measurements of pseudoscalar me¬ 
son photoproduction reactions with unprecedented qual¬ 
ity at facilities like ELSA, MAMI, and JLab have opened 
a path towards a more complete picture of the baryon 
spectrum. The photoproduction of r] mesons is isospin 
selective and allows for an analysis of N* states unaf¬ 
fected by contributions from A* states. Furthermore, the 
t]N final state is physically open for all resonances in the 
second resonance region and beyond. Eta photoproduc¬ 
tion is, thus, a prime reaction for resonance analysis and 
future complete experiments. Recently, polarization ob¬ 
servables with large angular coverage and high statistics 
have emerged. Among them are the target asymmetry T 
and the beam-target asymmetry F. The latter observ¬ 
able has been measured at MAMI for the very first time 
in 7 p rjp. 

However, even with the measurement of more observ¬ 
ables and an improved coverage of the data in angles 
and energy, a reliable determination of resonance prop¬ 
erties requires a combined analysis of reactions with dif¬ 
ferent initial and final states. One of these reactions is 
TT~p — >■ pn where, however, the data situation is known 
to be problematic. To avoid bias, this requires a refit not 
only of parameters tied to photoproduction, but also of 
hadronic parameters. Baryon resonance analyses would 
greatly benefit from a re-measurement of the tt~p —>■ pn 
reaction. 

Dynamical coupled-channel (DCC) approaches provide 
an especially suited tool to combine different reaction 
channels in a global analysis. In the present study, we 
extended the Jiilich DCC framework to eta photoproduc¬ 
tion. Based on a simultaneous analysis of nearly 30,000 
data points for pion and eta photoproduction off the pro¬ 
ton and the world database on the pion-induced reactions 
ttN —>■ ttN, pN, KA, and 77E, we extracted the spec¬ 
trum of nucleon and A resonances in terms of pole posi¬ 
tions, residues and photocouplings at the pole in an en- 
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ergy regime from ttN threshold up to E ~2.3 GeV. In the 
current approach, unitarity and analyticity are respected 
which is a prerequisite for a reliable determination of the 
resonance properties. 

Poles and residues were compared to the preceding 
Julich2012 analysis [10] in which only hadronic data 
were considered. The effect of the photoproduction data 
is most apparent for higher resonances, but also no¬ 
ticeable in case of well established states like the two 
S'!! resonances iV(1535) 1/2“ and A^(1650) 1/2“ whose 
widths change when photoproduction data are included. 
Also, some photocouplings at the pole changed in the 
present analysis compared to the Julich2013 solution [51] 
in which only pion photoproduction data, but no eta pho¬ 
toproduction data were considered. 

In order to estimate the influence of the recent MAMI 
T and F measurements, two different fits were performed, 
including the new data only in the second fit. Changes 
in the resonance parameters are predominantly observed 
for less well established states like the A^(1710) 1/2+ or 
higher lying resonances. Smaller but significant changes 
appear also for well reputed states and particularly for 
the photocouplings at the pole. Moreover, the new data 
on T and F have a major influence on the multipoles. 

In general, the multipole content of eta photoproduc¬ 
tion is less well established than for pion photoproduc¬ 
tion. This calls for further measurements of single- and 
double polarization observables. Upcoming experiments 
on polarization observables will have significant impact 
on the resonance spectrum and will help to identify so- 
called missing states and determine their resonance pa¬ 
rameters. 
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Appendix A: Renormalization of the nucleon mass 
and coupling 

The nucleon is included as an s-channel state in the 
Pii partial wave. In contrast to the other explicit states 
in this partial wave the bare mass and the bare cou¬ 
pling of the nucleon are not free parameters but un¬ 
dergo a renormalization process such that the nucleon 
pole position and residue to the nN channel correspond 


to the physical values, i.e. Eq = = 938 MeV and 

firNN = 0.964 [95]. Note that in the present study the 
nucleon is only allowed to couple to the ttN channel. Ef¬ 
fects of the coupling to other channels with significantly 
higher threshold energies are small and can be absorbed 
in the renormalization process. 

In Ref. [10] the renormalization of the nucleon in the 
presence of two s-channel states was illustrated. In the 
present study, we introduce an additional contact term 
in the Pn partial wave. Hence, the renormalization pro¬ 
cedure has to be modified. 

For this purpose we define the following reduced self¬ 
energies E 

_ / pb \2^red _ 2 p2 v^red. _ 2 

11 — Un) ^11 — ^ JttNN ^11 — ^ ^11 
Ei 2 = 

S2I = =xf^NNT.-^21 =X^21 

ei3 = = = 

E31 = fl S-d = X Umn = a; E31 , (Al) 

with a; € K and the bare ttNN coupling constant /^, 

In = X fnNN ■ (A2) 


We also define the reduced nucleon resonance vertices 




via 


-piOjC 

i /i;l 


rb -pred;a,c 
JN^ 


X JttNN P 


red;a 


,c 


— rj. f 


(A3) 


The nucleon pole position at Eq = is given by 

a zero of the determinant of D~^ of Eq. (10). Using 
Eqs. (Al) and (A3) we obtain an expression for the bare 
nucleon mass 


_ ™Phys 


rrijy = m 


N 


+ X^ 


- E 


11 


(A4) 


-f 


G-^E^2 + Si 3 (G^^Ei 3 + 2E12E23) 


(Ei 3 - G-^G3“ 


) 


Here, we introduced the auxiliary quantities G 2 ^ and 
G 3 -': 

G2-1 = G^\E) ={E-m\- E22(E)) 

G3“1 = G3“i(U) = m^-E33(E)|^^^. . (A5) 

phys 

In Eq. (A4) all quantities E(U), G^^(U), and G'^^{E) 
are evaluated at U = 

To determine x and thus the bare nucleon coupling /^, 
we exploit that at the nucleon pole the physical residue 
{a-i)TrN^TrN has to agree with the residue of T+ from 
Eq. (9). The physical residue is given by 


— Tl Tl : (-^^) 

where 7 / ( 7 /) are the bare nucleon vertices calculated at 
E = Eq with the physical nucleon coupling /,riv n instead 
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of the bare coupling (cf. Appendix B.l. of Ref. [50]), with Be ■= -§e- We obtain 


7i = 


3 f-^NN k En + + tun 

8 TT TTItt E]\fUJTr{EN + TTIn) 


(A7) 


where k is the particle momentum in the center-of-mass 
frame. The residue of of Eq. (9) can be calculated as 


R,eSp _™physT' 

h/Q—m pj 


detH-i 

dEdetD~^ 


E—m 


phys 

N 


(A 8 ) 

I 


detH 


7i 7i = 


Se det D 




E — m\ — 

-xEi2 

-xEi3 \ ' 




E -m\- S22 

-S23 1 

i ■ 

(A9) 

-XE31 

— E32 

ttim — E33/ 

Um;3/ 



Both sides of Eq. (A9) are evaluated at if = . Solv¬ 

ing Eq. (A9) we arrive at an expression for x which de¬ 
pends only on known or fitted quantities. We can calcu¬ 
late the bare mass and coupling of the nucleon by insert¬ 
ing this expression for x in Eqs. (A2) and (A4). Setting 


E 33 = Ei 3 = S 23 = 0 we recover the two resonance case, 
cf. Sec. 2.2 of Ref. [10]. 

The renormalization procedure is performed for each 
step in the fitting process. 
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TABLE II: Properties of the I = 1/2 resonances: Pole positions Eq (Ptot defined as -2Im_Bo), elastic ttN residues 
(I^TrivI, 07 rAr-> 7 rAr), and the normalized residues (yTVNT^/Ptot, of the inelastic reactions -kN /r with /r = tjA", KK, 

KT,. (*): not identified with PDG name; (a): dynamically generated. The resonance parameters of fit Ahad were extracted 
from fit A of Ref. [10], where only hadronic data were considered. 


fit 

Re Eo 

[MeV] 

-2Im Eo 

[MeV] 

[r.,riv[ 

[MeV] 

[deg] 

-pl/2pl/2 

^ 7tN^ r]N 

O'K N —^r]N 

[deg] 

pl/2pl/2 

^ KA 

d-irN^KA 

[deg] 

pl/2pl/2 

d-nN^KS 

[deg] 

Ptot 

]%] 

Ptot 

[%] 

Ptot 

[%] 

A(1535) 1/2“ 

A 

1497 

105 

23 

-48 

51 

110 

5.6 

-26 

8.9 

-78 


B 

1499 

104 

22 

-46 

51 

112 

5.0 

32 

5.0 

-69 


Ahad 

1498 

74 

17 

-37 

51 

120 

7.7 

68 

15 

-74 

A(1650) 1/2" 

A 

1664 

126 

31 

-59 

20 

47 

21 

-56 

26 

-76 


B 

1672 

137 

37 

-59 

21 

48 

20 

-54 

26 

-74 


Ahad 

1677 

146 

45 

-43 

15 

57 

25 

-46 

26 

-63 

iV(1440) 1/2+) 

A 

1349 

221 

32 

-69 

5.0 

11 

4.5 

-178 

1.1 

140 


B 

1355 

215 

62 

-98 

7.8 

-27 

16 

145 

2.7 

113 


Ahad 

1353 

212 

59 

-103 

2 

-40 

2 

156 

1 

67 

iV(1710) 1/2+ 

A 

1611 

140 

2.7 

-40 

6.1 

175 

3.7 

-49 

0.9 

-58 


B 

1651 

121 

3.2 

55 

16 

-180 

12 

-32 

0.4 

-43 


Ahad 

1637 

97 

4 

-30 

24 

130 

9.4 

-83 

3.9 

-136 

N{1750) l/2+_^) 

A 

1746 

317 

10 

-133 

0.7 

-85 

0.3 

127 

1.8 

-48 


B 

1747 

323 

14 

-144 

0.2 

138 

0.4 

86 

1.6 

-55 


Ahad 

1742 

318 

8 

161 

0.5 

-140 

0.8 

-170 

2.2 

4 

A(1720) 3/2+ 

A 

1711 

209 

5.3 

-45 

1.0 

132 

0.3 

-97 

2.2 

96 


B 

1710 

219 

4.2 

-47 

0.7 

106 

1.1 

-70 

0.2 

79 


Ahad 

1717 

208 

7 

-76 

1.2 

98 

3.1 

-89 

1.7 

64 

A(1520) 3/2” 

A 

1512 

107 

42 

-15 

4.6 

82 

7.8 

155 

8.9 

140 


B 

1512 

89 

37 

-6 

2.6 

95 

6.9 

158 

4.9 

-41 


Ahad 

1519 

110 

42 

-16 

3.5 

87 

5.8 

158 

0.8 

163 

iV(1675) 5/2- 

A 

1649 

129 

25 

-21 

6.1 

-43 

0.3 

90 

2.0 

-172 


B 

1646 

125 

24 

-22 

4.4 

-43 

0.1 

100 

3.1 

-175 


Ahad 

1650 

126 

24 

-19 

6.0 

-40 

0.3 

-93 

3.3 

-168 

iV(1680) 5/2+ 

A 

1668 

107 

37 

-23 

1.0 

130 

0.7 

-96 

0.1 

145 


B 

1669 

100 

34 

-19 

2.7 

136 

0.1 

90 

0.4 

148 


Ahad 

1666 

108 

36 

-24 

0.4 

-47 

0.2 

-99 

0.1 

141 

A(1990) 7/2+ 

A 

1791 

305 

5.9 

-93 

0.6 

-106 

1.6 

-131 

1.2 

21 


B 

1738 

188 

4.3 

-70 

1.3 

-82 

2.2 

-111 

0.5 

24 


Ahad 

1788 

282 

4 

-84 

0.4 

-99 

1.7 

-123 

0.8 

28 

A(2190) 7/2- 

A 

2081 

359 

45 

-27 

0.4 

137 

2.0 

-48 

1.2 

-59 


B 

2074 

327 

35 

-40 

1.6 

129 

0.5 

-51 

1.3 

-69 


Ahad 

2092 

363 

42 

-31 

0.1 

-28 

1.9 

-51 

1.3 

-63 

iV(2250) 9/2- 

A 

2165 

487 

21 

-67 

0.7 

-94 

0.6 

-105 

0.4 

-115 


B 

2062 

403 

8.2 

-64 

1.7 

-89 

0.6 

-101 

0.2 

70 


Ahad 

2141 

465 

17 

-67 

0.6 

-92 

1.1 

-103 

0.3 

-114 

N{2220) 9/2+ 

A 

2173 

611 

70 

-62 

0.3 

-103 

1.0 

60 

0.1 

-136 


B 

2171 

593 

62 

-59 

0.4 

-101 

0.7 

62 

0.9 

-128 


Ahad 

2196 

662 

87 

-67 

0.1 

63 

0.9 

53 

0.8 

-138 
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TABLE III: Properties of the 7 = 3/2 resonances: Pole positions Eq (Ptot defined as -2Im7?o), elastic -rN residues 
(|r,rjv|, 07 rAr-> 7 rAr), and the normalized residues (•^rTrArPjj/Ptot, of the inelastic reactions nN —> KT, and rN rA. 

(a): dynamically generated. The resonance parameters of fit Ahad were extracted from fit A of Ref. [10], where only hadronic 
data were considered. 


fit 

Pole position 

Re Eq 

[MeV] 

-2Iin Eq 

[MeV] 

rN Residue 

[rTTN^] 

]MeV] [deg] 

KT> chani 

pl/2pl/2 

lel 

[deg] 

ttA, chan 

pl/2pl/2 

-kN^ ttA 

nel (6) 

^7riV—^ ttA 

[deg] 

ttA, chan 

pl/2pl/2 

^ ttAT-^ ttA 

nel (7) 

^TT/V—>-77 A 

[deg] 

Ptot 

[%] 

Ptot 

[%] 

Ptot 

[%] 

A(1620) 1/2“ 

A 

1599 

69 

17 

-106 

21 

-106 

- 

- 

57 

103 


B 

1600 

65 

16 

-104 

22 

-105 

- 

- 

57 

105 


Ahad 

1599 

71 

17 

-107 

22 

-107 

- 

- 

57 

102 

A(1910) 1/2+ 

A 

1799 

651 

83 

-83 

2.1 

-129 

54 

130 

- 

- 


B 

1799 

648 

90 

-83 

1.9 

-123 

58 

131 

- 

- 


Ahad 

1788 

575 

56 

-140 

4.7 

-144 

41 

71 

- 

- 

A(1232) 3/2+ 

A 

1218 

91 

45 

-36 








B 

1218 

92 

46 

-36 








Ahad 

1220 

86 

44 

-35 







A(1600) 3/2+) 

A 

1552 

348 

24 

-156 

14 

-6.1 

33 

31 

1.4 

36 


B 

1552 

350 

23 

-155 

13 

-5.6 

31 

31 

1.3 

29 


Ahad 

1553 

352 

20 

-158 

11 

-7 

28 

28 

1 

15 

A(1920) 3/2+ 

A 

1714 

882 

36 

147 

15 

-34 

6.3 

132 

1.3 

-115 


B 

1715 

882 

38 

146 

17 

-35 

6.9 

131 

1.3 

-115 


Ahad 

1724 

863 

36 

163 

16 

-21 

7 

144 

1 

-101 

A(1700) 3/2“ 

A 

1676 

305 

24 

-8.4 

1.3 

-148 

5.4 

166 

39 

150 


B 

1677 

305 

24 

-7.3 

1.1 

-147 

5.4 

166 

39 

151 


Ahad 

1675 

303 

24 

-9 

1.5 

-150 

5 

166 

39 

149 

A(1930) 5/2" 

A 

1797 

655 

22 

-155 

3.3 

0.3 

14 

30 

0.6 

131 


B 

1836 

724 

34 

-155 

4.3 

-0.5 

15 

30 

0.9 

121 


Ahad 

1775 

646 

18 

-159 

3.1 

-3 

12 

26 

< 1 

- 

A(1905) 5/2+ 

A 

1795 

245 

5.3 

-82 

0.2 

-148 

0.8 

85 

9.7 

80 


B 

1795 

247 

5.3 

-89 

0.1 

-155 

0.9 

64 

8.7 

72 


Ahad 

1770 

259 

17 

-59 
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TABLE IV: Properties of the 7 = 1/2 (left) and I = 3/2 resonances (right): photocouplings at the pole (Ap^j^, i9^) according 
to Eq. (15). (*): not identified with PDG name; (a): dynamically generated. The photocouplings labeled “fit 2” were extracted 
in Ref. [51], where parameters of the hadronic T matrix were not altered, i.e. the photocouplings of fit 2 correspond to the pole 
positions and residues of fit Ahad in Tables II and III. 
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